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Abstract. We show that under the definition of computabihty which is natural from 
the point of view of applications, there exist non-computable quadratic Julia sets. 
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. . 1. Summary of the paper 

^. 

r*\ . Polynomial Julia sets have emerged as the most studied examples of fractal sets generated 

by a dynamical system. Apart from the beautiful mathematics, one of the reasons for their 

popularity is the beauty of the computer-generated images of such sets. The algorithms 

used to draw these pictures vary; the most naive work by iterating the center of a pixel 

to determine if it lies in the Julia set. Milnor's distance estimator algorithm |Milj . uses 

classical complex analysis to give a one-pixel estimate of the Julia set. This algorithm and 

^ ' its modifications work quite well for many examples, but it is well known that in some 

^ ■ particular cases computation time will grow very rapidly with increase of the resolution. 

Moreover, there are examples, even in the family of quadratic polynomials, when no satis- 

^ I factory pictures of the Julia set exist. In this paper we study computabihty properties of 

2 ■ Julia sets of quadratic polynomials. Under the definition we use, a set is computable, if, 

(^ . roughly speaking, its image can be generated by a computer with an arbitrary precision. 

"^ ! Under this notion of computabihty we show: 

c^ ■ Main Theorem There exists a parameter value c G C such that the Julia set of the 

quadratic polynomial fc{z) = z^ + c is not computable. 

The structure of the paper is as follows. In the Introduction we discuss the question of 
^ ■ computabihty of real sets, and make the relevant definitions. Further in this section we 

briefly introduce the reader to the main concepts of Complex Dynamics, and discuss the 
properties of Julia sets relevant to us. In the end of the Introduction, we outline the 
conceptual idea of the proof of Main Theorem. Section §3 contains the technical lemmas 
on which the argument is based. In §4 we complete the proof. 
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2. Introduction 
2.1. Introduction to the Computability of Real Sets. 

Classical Computability. The computability theory in general allows us to classify prob- 
lems into the tractable ("computable") and intractable ("uncomputable"). All the common 
computational tasks such as integer operations, list sorting, etc. are easily seen to be com- 
putable. On the other hand, there are many uncomputable problems. 

In the formal setting for the study of computability theory computations are performed 
by objects called the Turing Machines. Turing Machines were introduced in 1936 by Alan 
Turing (see |Tur j ) and are accepted by the scientific community as the standard model of 
computation. The Turing Machine (TM in short) is capable of solving exactly the same 
problems as an ordinary computer. Most of the time, one can think of the TM as a 
computer program written in any programming language. It is important to mention that 
there are only countably many TMs, which can be enumerated in a natural way. See |Sip| 
for a formal discussion on TMs. We define computability as follows. 

Definition 2.1. We say that a function / : {0, 1}* -^ {0, 1}* is computable if there is a 
TM, which on input string s outputs the string /(s). 

We say that the set L C {0, 1}* is computable or decidable if its characteristic function 
Xl '■ {0, 1}* — > {0, 1} is computable. 

While most "common" functions are computable, there are uncountably many uncom- 
putable functions and undecidable sets. The best-known intractable problems are the 
Halting Problem and the solvability of a Diophantine equation (Hilbert's 10-th problem), 
see |Sip| and |Matj for more information. 



Computability of Real Functions and Sets. In the present paper we are interested 
in the computability of functions / : M" ^ M and subsets of M", particularly subsets of 
M^ = C. We cannot directly apply Definition 12.11 here, since a real number cannot be 
represented in general by finite sequences of bits. 

Denote by D the set of the dyadic rationals, that is, rationals of the form ^. We say 
that : N ^ D is an oracle for a real number x, if \x — 4>{n)\ < 2~" for all n eN. In other 
words, (p provides a good dyadic approximation for x. We say that a TM M'^ is an oracle 
machine, if at any step of the computation M is allowed to query the value 0(n) for any 
n. This definition allows us to define the computability of real functions on compact sets. 

Definition 2.2. We say that a function / : [a,b] — > [c,d] is computable, if there exits an 
oracle TM M'^{rn) such that if (j) is an oracle for x G [a, 6], then on input m, M"^ outputs 
a 1/ e D such that \y - f{x)\ < T"^ . 

This definition was first introduced by Grzegorczyk JGrzj and Lacombe |Lacj . and follows 
in the tradition of Computable Analysis originated by Banach and Mazur in 1937 (see 

jMii). 
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To understand the definition better, the reader without a Computer Science background 
should think of a computer program with an instruction 

READ real number x WITH PRECISION n{m). 

On the execution of this command, a dyadic rational d is input from the keyboard. This 
number must not differ from x by more than 2""'^™'' (but otherwise can be arbitrary). The 
algorithm then outputs f{x) to precision 2"". 

In other words, with an access to arbitrarily good approximations for x, M should be able 
to produce an arbitrarily good approximation for f{x). This definition trivially generalizes 
to domains of higher dimension. See |Kolj for more details. One of the most important 
properties of computable functions is that 

Proposition 2.1. Computable functions are continuous. 

Let fi" C M'^ be a compact set. We would like to give a definition for K being computable. 
Note that saying by a naive analogy with Definition 12.11 that K is computable if and only 
if the characteristic function xk is computable does not work here, since by the above 
proposition only continuous functions can be computable. 

We say that a TM M computes the set K if it approximates K in the Hausdorff metric. 
Recall that the Hausdorff metric is a metric on compact subsets of R" defined by 

(2.1) dniX, Y) = inf {e > | X C U,{Y) and Y C U,{X)}, 

where U^{S) is defined as the union of the set of e-balls with centers in S. 
We introduce a class C of sets which is dense in metric dn among the compact sets and 
which has a natural correspondence to binary strings. Namely C is the set of finite unions 
of dyadic balls: 

r " 

C = < [J B{di, Ti) I where di, r^ G 1 

U=i 

The following definition is equivalent to the set computability definition given in |Weij (see 
also IMl). 

Definition 2.3. We say that a compact set i^ C M'^ is computable, if there exists a TM 
M{m), such that on an input m & N, the machine M{m) outputs an encoding of Cm G C 
such that dniK, Cm) < 1"^ ■ 

To illustrate the robustness of this definition we present the following two equivalent 
characterizations of computable sets. The first one relates the definition to computer 
graphics. It is not stated precisely here, but it can be easily made precise. The second one 
relates the computability of sets to the computability of functions as per Definition 12.21 

Theorem 2.2. For a compact K (Z^ the following are equivalent: 

(1) K is computable as per Definition \2.!A 

(2) (in the case k = 1,2) K can be drawn on a computer screen with arbitrarily good 
precision, 
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(3) the distance function dxix) = inf{|x — y\ \ y E K} is computable as per Definition 



In the present paper we are interested in questions concerning the computabihty of the 
Juha set Jc = J{fc) = J{z^ + c) (see the next section for the definition). Since there 
are uncountably many possible parameter values for c, and only countably many TMs, we 
cannot expect for each c to have a machine M such that M computes Jc- On the other 
hand, it is reasonable to want M to compute Jc with an oracle access to c. Define the 
function J : C ^ K* [K* is the set of all compact subsets of C) by J{c) = J{fc)- In a 
complete analogy to Definition 12.21 we can define 

Definition 2.4. We say that a function f : S ^ K* for some bounded set S C M'^ is 
computable, if there exits an oracle TM M'^{m) such that if (p is an oracle for x & S, then 
on input m, M'^ outputs a set Cm G C such that duiCm, f{x)) < 2""*. 

In the case of Julia sets: 

Definition 2.5. We say that Jc is computable if the function J : d ^-^ Ja is computable 
on the set {c}. 

The following has been shown (see jBrvlj . |Ret j ) : 

Theorem 2.3. Denote by Ti the set of parameters c for which Jc is hyperbolic, then 
(i) Jc is computable for all c ^Ti, moreover 
(a) the function J is computable on each bounded subset ofTi. 

Our goal in this paper is to show that there are values of c for which Jc is not computable 
under Definition 12.51 which is the weakest possible definition in this setting. We will be 
using the following version of Theorem 12.11 for set functions. 

Theorem 2.4. Suppose that a TM M^ computes the function J on a set S . Then J is 
continuous on S in Hausdorff sense. 

Proof. Let c be any point in S, and let e = 2~^ be given. Let be an oracle for c such 
that \(j){n) — c\ < 2^("+^) for all k. We run M'^{k + 1) with this oracle 0. By the definition 
of J, it outputs a set L which is a 2^^^+^-' approximation of Jc in the Hausdorff metric. 

The computation is performed in a finite amount of time. Hence there is an m such that 
is only queried with parameters not exceeding m. Then for any x such that \x—c\ < 2~(™-+^), 
is a valid oracle for x up to parameter value of m. In particular, we can create an oracle 
■0 for X that agrees with on 1, 2, ... ,771. If x G 5, then the execution of M^{k + 1) 
will be identical to the execution of M'^{k + 1), and it will output L which has to be an 
approximation of J^c- Thus we have 

dHiJc, J.) < dniJc, L) + rf^(4, L) < 2-^^+^^ + 2-^^+^^ = 2"^ 

This is true for any x & B{c, 2"^™+^^) fl S. Hence J is continuous on S. D 
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In the next section we proceed to define Julia sets of rational maps and review their basic 
properties. In particular, towards the end of the introduction, we will see a mechanism by 
which the continuity required by Theorem 12.41 mav fail. 

It should be noted that the question of computability of dynamically generated fractal 
sets, such as Julia sets, has been discussed by Blum, Cucker, Shub, and Smale in jBCSSj . 
The definition of set computability used in |BCSSj is, however, quite different from Defi- 
nition 12.31 The BCSS model allows infinite-precision arithmetic, but requires completely 
accurate pictures to be generated. Under this definition all Julia sets but the most trivial 
ones can be shown to be uncomputable. 

2.2. Julia sets of polynomial mappings. We recall the main definitions of complex 
dynamics relevant to our result only briefly; a good general reference is the book of Milnor 
|Milj . For a rational mapping R of degree deg R = d > 2 considered as a dynamical system 
on the Riemann sphere 

R:C-^C 
the Julia set is defined as the complement of the set where the dynamics is Lyapunov-stable: 

Definition 2.6. Denote F{R) the set of points 2; G C having an open neighborhood U{z) 
on which the family of iterates R"'\u{z) is equicontinuous. The set F{R) is called the Fatou 
set of R and its complement J{R) = C \ F{R) is the Julia set. 

In the case when the rational mapping is a polynomial 

P{z) = ao + aiz + ■ ■ ■ + adz'^ : C ^ C 

an equivalent way of defining the Julia set is as follows. Obviously, there exists a neigh- 
borhood of cxD on C on which the iterates of P uniformly converge to cxd. Denoting A{oo) 
the maximal such domain of attraction of cxd we have A{oo) C F{R)- We then have 

J{P) = dA{oo). 

The bounded set C\cly4(oo) is called the filled Julia set, and denoted K{P); it consists of 
points whose orbits under P remain bounded: 

K{P) = {zeC\ sup |P"(2)| < cx)}. 

n 

For future reference, let us list in a proposition below the main properties of Julia sets: 

Proposition 2.5. Let R : C ^ C be a rational function. Then the following properties 
hold: 

• J{R) is a non-empty compact subset ofC which is completely invariant: R'^{J{R)) = 
J{R); 

• J{R) = J{R") for all n E N; 

• J{R) is perfect; 

• if J{R) has non-empty interior, then it is the whole ofC; 
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• let U C C be any open set with U fl J{R) ^ 0. Then there exists n E N such that 
R'iU) D J{R); 

• periodic orbits of R are dense in J{R). 

Let us further conunent on the last property. For a periodic point zq = R^{zq) of period 
p its multiplier is the quantity A = A(zo) = DRp{zq). We may speak of the muhipher of a 
periodic cycle, as it is the same for all points in the cycle by the Chain Rule. In the case 
when I A| 7^ 1, the dynamics in a sufficiently small neighborhood of the cycle is governed by 
the Mean Value Theorem: when |A| < 1, the cycle is attracting {super- attracting if A = 0), 
if |A| > 1 it is repelling. Both in the attracting and repelling cases, the dynamics can be 
locally linearized: 

(2.2) ,l,{R^{z)) = \-^^{z) 

where ^^ is a conformal mapping of a small neighborhood of zq to a disk around 0. By 
a classical result of Fatou, a rational mapping has at most finitely many non-repelling 
periodic orbits. Therefore, we may refine the last statement of Proposition 12.51 

• repelling periodic orbits are dense in J{R). 

In the case when |A| = 1, so that A = e^'^*^, ^ G M, the simplest to study is the parabolic 
case when 6 = n/m G Q, so A is a root of unity. In this case i?^ is not locally linearizable; 
it is not hard to see that zq £ J{R). In the complementary situation, two non- vacuous 
possibilities are considered: Cremer case, when R^ is not linearizable, and Siegel case, 
when it is. In the latter case, the linearizing map if) from (j2.2j) conjugates the dynamics 
of R^ on a neighborhood U{zq) to the irrational rotation by angle 9 (the rotation angle) 
on a disk around the origin. The maximal such neighborhood of zq is called a Siegel disk. 
Siegel disks will prove crucial to our study, and will be discussed in more detail in the next 
section. 

To conclude the discussion of the basic properties of Julia sets, let us consider the sim- 
plest examples of non-linear rational endomorphisms of the Riemann sphere, the quadratic 
polynomials. Every affine conjugacy class of quadratic polynomials has a unique represen- 
tative of the form fc{z) = z^ + c, the family 

f,{z) = z^ + c, ceC 

is often referred to as the quadratic family. For a quadratic map the structure of the Julia 
set is governed by the behavior of the orbit of the only finite critical point 0. In particular, 
the following dichotomy holds: 

Proposition 2.6. Let K = K{f^ denote the filled Julia set of fc, and J = J{fc) = dK . 
Then: 

• G -ft' implies that K is a connected, compact subset of the plane with connected 
complement; 

• ^ K implies that K = J is a planar Cantor set. 

The Mandelbrot set A4 G C is defined as the set of parameter values c for which J{fc) is 
connected. 
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Continuity of the dependence c i-^ Jifc)- A natural question to pose for polynomials 
in the quadratic family is whether the Julia set varies continuously with the parameter 
c. To make sense of this question, recall the definition of the Hausdorff distance dist^:/ 
between compact sets X, Y in the plane ()2.H) . It turns out that the dependence c i— > J(/c) 
is discontinuous in the Hausdorff distance. For an excellent survey of this problem see 
the paper of Douady |Doj . The discontinuity which has found most interesting dynamical 
applications occurs at parameter values for which fc has a parabolic point. We, however, 
will employ a more obvious discontinuity which is related to Siegel disks. Let us first note 
that by a result of Douady and Hubbard ;DHlj a quadratic polynomial has at most one 
non-repelling cycle in C. In particular, there is at most one cycle of Siegel disks. 

Proposition 2.7. Let c^ & Ai be a parameter value for which fc has a Siegel disk. Then 
the map c i-^ J{fc) is discontinuous at c*. 

Proof. Let zq be a Siegel periodic point of fc and denote A the Siegel disk around Co, V 
its period, and 9 the rotation angle. By the Implicit Function Theorem, there exists a 
holomorphic mapping C, : f/(c*) -^ C such that C(c*) = Zq and C(c) is fixed under {fcY- 
The mapping 

v.c^ DifcnCic)) 

is holomorphic, hence it is either constant or open. If it is constant, all quadratic polyno- 
mials have a Siegel disk. This is not possible: for instance, /1/4 has a parabolic fixed point, 
and thus no other non-repelling cycles. Therefore, u is open, and in particular, there is 
a sequence of parameters c„ — »• c* such that C(cn) has multiplier e^'^*^"/^". Since C(cn) is 
parabolic, it lies in the Julia set of /c„ . Hence 

distH(J(/,J, J(/,J) > dist(c„9A)/2 

for n large enough. D 

Thus an arbitrarily small change of the multiplier of the Siegel point may lead to an 
implosion of the Siegel disk - its inner radius collapses to zero. We make a note of an 
immediate consequence of the above proposition and Theorem 12.41 

Proposition 2.8. For any TM M'^{n) with an oracle for c G C denote Sm the set of all 
values of c for which M'^ computes Jc- Then Sm 7^ C. 

In other words, a single algorithm for computing all quadratic Julia sets does not exist. 



Siegel disks of quadratic maps. Let us discuss in more detail the occurrence of Siegel 
disks in the quadratic family. For a number G [0,1) denote [ri,r2, . . . ,r„, . . .], r^ G 



8 M. BRAVERMAN, M. YAMPOLSKY 

N U {00} its possibly finite continued fraction expansion: 
(2.3) [ri,r2,...,r„,...] = 



1 

ri H 



1 

r2 H 

••• + 



Sucli an expansion is defined uniquely if and only if ^ ^ Q. In this case, the rational 
convergents Pn/Qn = [1^1, ■ ■ ■ ,^n] are the closest rational approximants of 6 among the 
numbers with denominators not exceeding g„. In fact, setting A = e^'^*^, we have 

|A^ - 1| > lA*^" - 1| for alio < h < g„+i, h ^ g„. 

The difference |A''" — 1| lies between 2/qn+i and 27r/g„+i, therefore the rate of growth of 
the denominators g„ describes how well 9 may be approximated with rationals. 

Definition 2.7. The diophantine numbers of order k, denoted I^(fc) is the following class 
of irrationals "badly" approximated by rationals. By definition, 6 G T>{k) if there exists 
c > such that 

Qn+i < cq'^~^ 

The numbers g„ can be calculated from the recurrent relation 

g„+i = rn+iqn + qn-i, with go = 0, gi = 1. 

Therefore, 6 G 'D(2) if and only if the sequence {tj} is bounded. Dynamicists call such 
numbers bounded type (number-theorists prefer constant type). An extreme example of a 
number of bounded type is the golden mean 

^, = ^^ = [1,1,1,...]. 

The set 

V{2+) = f]V, 

fc>2 

has full measure in the interval [0, 1). In 1942 Siegel showed: 



Theorem 2.9 ( |Siej ). Let R be an analytic map with an periodic point zq E C of period p. 
Suppose the multiplier of the cycle 

A = e^^*'' with e e V{2+), 



then the local linearization equation h2.2\) holds. 

The strongest known generalization of this result was proved by Brjuno in 1972: 
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Theorem 2.10 ( |Bruj ). Suppose 

log(gn+i) 



(2.4) B{e) = J2 



< oo, 



Qn 

then the conclusion of Siegel's Theorem holds. 

Note that a quadratic polynomial with a fixed Sigel disk with rotation angle 6 after an 
affine change of coordinates can be written as 

(2.5) Pe{z) =z^ + e^^^'z. 

In 1987 Yoccoz |Yocj proved the following converse to Brjuno's Theorem: 

Theorem 2.11 ( |Yocj ). Suppose that for 9 G [0, 1) the polynomial Pq has a Siegel point at 
the origin. Then B{6) < oo. 

The numbers satisfying (|2.4j) are called Brjuno numbers; the set of all Brjuno numbers will 
be denoted B. It is evident that UV{k) C B. The sum of the series (|2.4p is called the 
Brjuno function. For us a different characterization of B will be more useful. Inductively 
define 6\ = and 6n+i = {l/6'„}. In this way, 

We define the Yoccoz' s Brjuno function as 



$(^) = V^1^2---^n-ll0g^. 

One can verify that 

B{e) < oo^<l>(0) < oo. 

The value of the function $ is related to the size of the Siegel disk in the following way. 

Definition 2.8. Let P{9) be a quadratic polynomial with a Siegel disk A^ 9 0. Consider 
a conformal isomorphism : © h-^ A fixing 0. The conformal radius of the Siegel disk Ag 
is the quantity 

r{9) = Wm. 
For all other 9 G [0, oo) we set r{9) = 0. 

By the Koebe One-Quarter Theorem of classical complex analysis, the internal radius of 
Ag is at least r{9)/4. Yoccoz |Yocj has shown that the sum 

^9) + \ogr{9) 

is bounded from below independently of 6* G i3. Recently, Buff and Cheritat have greatly 
improved this result by showing that: 

Theorem 2.12 ( |BC2j ). The function 9 i— >• $(6^) + log r(6') extends to R as a 1-periodic 
continuous function. 
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We remark that the following stronger conjecture exists (see |MM Yj ) : 

Marmi-Moussa-Yoccoz Conjecture. (MMYj ThefunctionO i— >• $(^)+logr(^) is Holder 
of exponent 1/2. 

Dependence of the conformal radius of a Siegel disk on the parameter. In this 
section we will show that the conformal radius of a Siegel disk varies continuously with the 
Julia set. To that end we will need a preliminary definition: 

Definition 2.9. Let (f/„, m„) be a sequence of topological disks Un C C with marked points 
Un G Un- The kernel OT Caratheodory convergence {Un,Un) — ^ {U,u) means the following: 

• for any compact K <ZU and for all n sufficiently large, K G Un', 

• for any open connected set W 3 u, ii W G Un for infinitely many n, then W G U. 

The topology on the set of pointed domains which corresponds to the above definition of 
convergence is again called kernel or Caratheodory topology. The meaning of this topology 
is as follows. For a pointed domain (f/, u) denote 

the unique conformal isomorphism with 0{c/,„)(O) = u, and (0(c/,u))'(O) > 0. We again 
denote r{U,u) = |(0(c/,n))'(O)| the conformal radius of U with respect to u. 
By the Riemann Mapping Theorem, the correspondence 

L : {U, u) ^ (f)(^u,u) 

establishes a bijection between marked topological disks properly contained in C and uni- 
valent maps : D — > C with 0'(O) > 0- The following theorem is due to Caratheodory, a 
proof may be found in |Pomj : 

Theorem 2.13 (Caratheodory Kernel Theorem). The mapping l is a homeomorphism 
with respect to the Caratheodory topology on domains and the compact-open topology on 
maps. 

Proposition 2.14. The conformal radius of a quadratic Siegel disk varies continuously 
with respect to the Hausdorff distance on Julia sets. 

Proof. To fix the ideas, consider the family Pe with 9 E B and denote A51 the Siegel 
disk of Pq. It is easy to see that the Hausdorff convergence J{Pe„) -^ J{.Pe) implies the 
Caratheodory convergence of the pointed domains 

(A,„,0)^(A,0). 

The proposition follows from this and the Caratheodory Kernel Theorem. D 

In fact, we can state the following quantitative version of the above result. For the proof, 
based on Koebe Theorem, see e.g. |RZj : 
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Lemma 2.15. Let U be a simply- connected bounded subdomain ofC containing the point 
in the interior. Suppose V C U is a simply- connected subdomain ofU, and dV C U^{dU). 
Then 

< r(f/, 0) - r(y, 0) < 4v/r(f/,0)v^. 

An immediate corollary is: 

Corollary 2.16. Suppose the function r{6) is uncomputable on the set {60} ■ Then the 
function 9 h-> J{Pe) is also uncomputable at the same point. 

Proof. Assume that J{Peo) is computable. Using the output of the TM computing this 
Julia set in an obvious way, for each e > we can obtain a domain V (zC such that 

V C Ae„ and dnidV, dAg^) < e. 

By Schwarz Lemma, the conformal radius r(^o) < 2. Hence, by Lemma f2. 151 

\r{V,0)-r{9o)\ < 6 = 8^^. 

Using any constructive version of the Riemann Mapping Theorem (see e.g. JBBJ), we can 
compute r{V, 0) to precision 6, and hence know r{6o) up to an error of 26. Given that 5 
can be made arbitrarily small, we have shown that r(^o) is computable. 

D 

We also state for future reference the following proposition: 

Proposition 2.17. Let {6i} be a sequence of Brjuno numbers such that 6i —>■ 9 and 
lim r{6i) = I > 0. Then 9 is also a Brjuno number and r{9) > I. 

Proof. Denote (pi = 0(Ae-,o)- Note that by Schwarz Lemma, the inverse ipi = {(pi)"^ linearizes 
Pg- on Ag- . By passing to a subsequence we can assure that 0j — > (p locally uniformly, and 
0'(O) > /. By continuity, (p~^ is a linearizing coordinate for Pg, so ^ is a Brjuno number. 
Moreover, 0(©) C Ag, and so by Schwarz Lemma r{9) > I. D 

Non-computability of the Yoccoz's Brjuno function. In addition to the non-compu- 
tability of the conformal radius, we also prove a non-computability result for the Yoccoz's 
Brjuno function $: 

Theorem 2.18 (Non-computability of $). There exists a parameter value 9 G R/Z 
such that $(0) < 00, and $(6') is not computable by any Turing Machine with an oracle 
for 9. 

It is worth noting that Marmi-Moussa-Yoccoz Conjecture as stated above and Theorem 
12.181 imply the existence of a non-computable quadratic Julia set. To see this, we first 
formulate: 

Conditional Implication 2.19. // Marmi-Moussa-Yoccoz Conjecture holds, then the 
function 

v. 9^(l>{9)+\ogr{9) 

is computable by one Turing Machine on the entire interval [0, 1]. 
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We use the following result of Buff and Cheritat ( |B(y2j ). 

Lemma 2.20 f |BC2p . For any rational point 6* = - G [0, 1] denote, as before, 

Pg{z) = e^^^'z + z\ 
and let the Taylor expansion of Pg''{z) at start with 

P°i(z) = z + Az'^+^ + . . . , forqeN 

Let L{6) = { — ] . Denote by ^tmnc the modification of $ applied to rational numbers 

where the sum is truncated before the infinite term. Then we have the following explicit 
formula for computing v{6): 

(2.6) v{e) = ^truuM + \ogL{e) + ^^^. 

Q 

Equation ()2.6|) allows us to compute the value oiv easily at every rational 9 G Qn[0, 1] with 
an arbitrarily good precision. In addition, assuming the conjecture, we have \v{x) —v{y) \ < 
2~" whenever |x — |/| < c ■ 2~^" for some constant c, hence v has an (easily) computable 
modulus of continuity. These two facts together imply that v is computable by a single 
machine of the interval [0, 1] (see for example Proposition 2.6 in |Ko2j ) . This implies the 
Conditional Implication. 
The following conditional result follows: 

Lemma 2.21 (Conditional). Suppose the Conditional Implication holds. Let 9 G [0,1] 
be such that ^{9) is finite. Then there is an oracle Turing Machine M{ computing $(^) 
with an oracle access to 9 if and only if there is an oracle Turing Machine Ml" computing 
r{9) with an oracle access to 9. 

Proof. Suppose that M{ computes ^{9) for some 9. Let M'^ be the machine uniformly 
computing the function v. Then we can use Mf and M'^ to compute logr(^) = v{9) — ^{9) 
with an arbitrarily good precision. We can then use this construction to give a machine 
Mg which computes r(9). 

The opposite direction is proved analogously. D 

Lemma [2.211 with Theorem 12. 181 implv that there is a for which r{9) is non-computable. 
Corollary 12.161 implies that for this value of 9 the Julia set of Pg is also non-computable. 
Note that for the proof of Conditional Implication we did not need the full power of the 
conjecture. All we needed is some computable bound on the modulus of continuity of v. 

Outline of the construction of a non-computable quadratic Julia set. We are now 

prepared to outline the idea of our construction. The outline given below is rather rough 
and suffers from obvious logical deficiencies. However, it captures the idea of the proof in a 
simple to understand form. Suppose that every Julia set of a polynomial Pg is computable 
by an oracle machine M'f', where represents 9. There are countably many machines, so 



NON-COMPUTABLE JULIA SETS 13 

we can enumerate them Mf,M2, . . .. Denote by Si the domain on which Mf computes 
Jpg properly. Then we must have: 

(1) c = u^i s., 

(2) for each i, the function J : 9 \-^ J{Pe) is continuous on Si. 




Figure 1. The Siegel disks of Pg for 6 given by the continued fractions 

[1, 1, 1, . . .], [1, 1, 1, 20, 1, . . .], and [1, 1, 1, 20, 1, 1, 1, 30, 1, . . .] 

Let us start with a machine M^ which computes J{PeJ) for 9^. = [1, 1, 1, . . .]. If any of 
the digits r^ in this infinite continued fraction is changed to a sufficiently large N E N, 
the conformal radius of the Siegel disk will become small. For N —>■ oo the Siegel disk 
will implode and its center will become a parabolic fixed point in the Julia set. Given the 
continuity of the dependence of the conformal radius of the Siegel disk on the Julia set, we 
have the following: 

There exists zi > 1 such that for every 9i whose continued fraction starts with zi ones, for 
the Julia set of Pq^ to be computable by M^^ , it must possess a Siegel disk of a conformal 
radius r(^i) > r(^,)(l - 1/8). 

We can thus "fool" the machine M^ by selecting 9i given by a continued fraction where 
all digits are ones except Tj^ = A^i >> 1. If we are careful, we can do it so that 

(2.7) r{9,)il - 1/4) < r{9i) < r{9,){l - 1/8). 

To "fool" the machine M^^ we then change a digit Vi^ for ^2 > H sufficiently far in the 
continued fraction of 9i to a large N2. In this way, we will obtain a Brjuno number ^2 for 
which 

(2.8) r(^,)(l - 1/4 - 1/8) < r(^2) < r(0*)(l - 1/4). 
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Continuing in this manner we will arrive at a limiting Brjuno number O^o for which the 
Julia set is uncomputable. To make such a scheme work, we need a careful analysis of 
the dependence of the conformal radius on the parameter. In this a key role is played by 
Theorem 12.121 of Buff and Cheritat which allows us to obtain a controlled change in the 
value of r{a) by changing $(«). The relevant analysis is carried out in the next section. 

Main analytic result. We formalize the strategy outlined above as follows: 

Theorem 2.22. There does not exist a partition of the circle R/Z into a countable union 
of sets Si such that for every i the function r{6) restricted to Si is continuous. 

The above formulation was suggested to us by John Milnor. Let us show how the Main 
Theorem and Theorem 12.181 follow from Theorem 12.221 

Proof of Main Theorem, assuming Theorem \2.2^ First we observe that there exists a pa- 
rameter 6*0 G M/Z such that the function r{9) is uncomputable on {^o}- Indeed, assume 
the contrary. There are only countably many Turing Machines with an oracle for 6. We 
enumerate them Mj, i G N in some arbitrary way (for instance, using the lexicographic 
order). Denote 

Si = {e e R/Z such that TM Mf computes r{e)} 
By Proposition 12. II the function r{6) is continuous on each of the Sj's, and we arrive at a 
contradiction with Theorem 12.221 

Now let us prove Main Theorem, again arguing by contradiction. Assume that for every 
c G C there exists a TM M*^ with an oracle for c which computes Jc- Let Pg = z"^ + e^'^^^ z 
as before. The affine change of coordinates transforming it into an element of the family 
fc is computable explicitly, and we have 

c = c{e) = AV4 - A/2 where A = e^'^^^ 

This implies that we can simulate an oracle for c given an oracle for 6. 

Set Co = c{9o) and consider the oracle TM M'f' computing the Julia set of /eg. By the 
above considerations, there exists an oracle TM M"^ with an oracle for 6 G R/Z which 
computes J{Pgg). This contradicts CoroUarv 12 . 1 61 and the proof is complete. 

D 

Proof of Theorem \2.1^ assuming Theorem \2.2'A Assume the contrary. Again, order in a 
sequence M/ , i G N all TMs with an oracle for 6 G R/Z. Let 

Vli = {6 e R/Z such that Uf computes the value of $(^)}. 

Denote VLq the set of all 9 with $(^) = oo. The value of r{d) on VLq is thus identically 0. 
Denote 

t;(^) = $(^)+logr(^), 

which by |BC2j continuously extends to R/Z. Given Proposition 12. H the function $(6') is 
continuous on each fij, z G N, and hence so is 

r{e) = exp{v{e) - ^{9)). 
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By our assumption, 

and we arrive at a contradiction with Theorem 12.221 

n 

A note on the connection with |B(ylj . A. Cheritat has pointed to us that methods of 
|BClj . where Siegel disks with smooth boundaries are constructed for the quadratic family 
can be used to derive the Main Theorem. We discuss this in the section following the 
proofs of the main theorems. We note here that the argument we give is based on quite 
elementary estimates of the function $ and is thus accessible to non-dynamicists. It has 
an added advantage of yielding Theorem 12.181 

3. Making Small Changes to $ 

3.1. Small Changes to $. A key step of the construction outlined above is making careful 
adjustments of r(6'j) as in the first two steps ()2.7p and ()2.8|) above. We do not have a direct 
control over the value of r{a), but Buff and Cheritat's Theorem 12.121 shows that small 
decreases of r{a) we would like to make correspond to a small controlled increment of the 
value of $(a). Estimates of a similar nature has appeared in the works of various authors 
(compare, for example, with |B(nj ). 
For a number 7 = [ai, 02, . . .] G M \ Q we denote 

«i(7) = 1 , 

ai H -, 



cti+i H ] 

ai+2 + ■ ■ ■ 

so that 

*(^) = X^ "i(7)"2(7) • • • ttn-i(7) log ^^ 

n>l 

The main goal of this section is to prove the following two lemmas: 



«n(7)' 



Lemma 3.1. For any initial segment I = [ai, 02, ... , a„], writeu = [ai, 02, ... , a„, 1, 1, 1, . . . 
Then for any e > 0, there is an m > and an integer N such that if we write (3^ = 
[oi, 02, ... , an, 1, 1, . . . , 1, A, 1, 1, . . .], where the N is located in the n + m-th position, and 

^{u) +£< <l>(/3^) < $(cj) + 2e. 

Lemma 3.2. For uo as above, for any e > there is an rriQ > 0, such that for any m > rriQ, 
and for any tail T = [an+m, Cin+m+i, ■ ■ ■] if we denote 

P = [Cti, 02, . . . , a„, 1, 1, . . . , 1, an+mi O-n+m+li ■■■]■, 

then 

$(/3^) > ^{uj) - e. 
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The proof is technical and will require some preparation. For Lemma \'A.1\ the idea is to 
choose an m large enough, so that changing a„+m (which will eventually be A^) by 1 changes 
the value of $ by a very small amount {< e). When N -^ oo, $(c<j) -^ oo, hence the value 
of $ must hit the interval {^{uj) + e, ^{co) + 2e). 
Denote 

$-(t^) = $(a;) - aia2 ■ ■ ■ an+m-i log . 

The value of the integer m > is yet to be determined. Denote 

P'^ = [ai, as, . . . , ttn, 1, 1, . . . , 1, iV, 1, 1, . . .]. 
We prove the following 
Lemma 3.3. For any N and i < n + m we have 



log 



ai{P 



N\ 



ai{(3 



N+l\ 



< 



oj— (n+m) 



N 



Proof. We prove the lemma by induction on i, starting from the base case i = n + m, and 
proceeding down to i = 0. The base case is i = n + m, we want to prove 



log 



a,+™(/3^+i) 



< 



A^' 



We have, 






where = (v^ + l) /2. Hence 



iV + 1/0 _A^ + 1 + 1; 



Ar+ 1 + 1/0 



iV+1/0 



iV+1/0' 



1 < ^n+m < 1 



N 



< e 



l/N 



and |logr„+m| < — . 

Induction step. Supposing that the statement is true for i + 1, we prove it for i. We have 

1 






Suppose that aj_|_i(/5^"'"^) > q;j+i(/3^). Then we know that 
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and we want to prove that 






since this expression is obviously bigger than 1. The situation is very similar in the case 
when aj+i(/?^"'"^) < q;j+i(/?^), with the numerator and the denominator exchanged. 

In other words, it is enough to prove that for < d < c < 1 and a pair of integers r > 0, 
k > 1 and a > 0, 

- < e" implies that < e"'^. 

a k + a 

First of all, it is easy to see that for k > 1, 



k+c 1+c 
k + d~ 1 + rf' 



hence it suffices to show that 



1 + c 



< e 



a/2 



(e") 



1/2 



< 



C\l/2 



1 + d 
Thus we need to demonstrate that 

1 + c 

1 + d ^ \dy 
This is equivalent to 

(1 + 2c + c'^)d < {l + 2d + d^)c ^ d + c^d<c + d^c ^ cd{c - d) < c - d. 

The last inequality holds, since cd < 1 and c — d > 0. 



n 



The following lemma is proven by induction exactly as the previous one with a different 
base. 

Lemma 3.4. Let 71 and 72 be two numbers whose continued fraction expansions coincide 
in the first n + m — 1 terms [ai, 02, ... , Qn+m-i] ■ Then for any i < n + m we have 



log 



a* (72 



< 



r)i— (?i+m)+l 



In particular, this applies with 71 = /5^ and 72 = P^- 

Proof. The proof goes by induction exactly as in Lemma 13.31 We need to verify the base 
case i = n + m — 1. For this value of i, 

Q;„+m-i(7i) = : , «„+„_! (72) = — 



with some /xi, /i2 £ [0, 1). Hence we have 

Oin+m-lill) 



'"n+m- 



-l+/i2 



log 



an+m-l{l2) 



< log2 < 1 = 2^ 



D 
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We now bound the influence of the difference on the log — terms. 



Oii 



Lemma 3.5. Let 'ji and '~f2 be two numbers whose continued fraction expansions coincide 
in the first n + m — 1 terms [oi, 02, . . . , a„+m-i] ■ Then for any i < n + m — 1 we have 

1 



log 



log- 



ai(7i) 



log 



1 



"1(72 



^ r)i— (n+m)+2 



In particular, this applies with 71 = /5^ and 72 = P^- 

Proof. Assume that 0:4(71) < 0:4(72), the reverse case is done in the same way. In this case 
we need to prove 

log^^^ 
\ — - < e 



log 



ai(72 



Denote c = Oj+i(7i) and d = 0^+1(72). Then we have Oj(7i) 
for some integer k > 1. Hence < d < c < 1. We have 

1 

log{k + c) 



and 0^(72) 



k + c 



k + d 



log 



log 



1 

ai(72 



log(/c + d) 



By Lemma Em we know that - < e 

d 



oi— {n+7Ti) + 2 



hence it suffices to show that 



log(A; + c) c 



< -. 



\og{k + d) d 

This is equivalent to < ; . Consider the function f(x) = on 

c d X 

the interval (0, 1). The reader can readily verify that f'{x) < for x G (0, 1) so that / is 

decreasing on this interval, and hence /(c) < f{d), which completes the proof. D 

We are now ready to bound the influence of changes in A^ on the value of <l>^. 

Lemma 3.6. For any uj of the form as in Lemma l,V. il and for any £ > 0, there is an 
rriQ > such that for any N and any m > mo, 

|$-(/?^)-<|.-(/3i)|<|. 

Proof. The ^ in the expression for $(/?^) converges, hence there is an mi > 1 such that the 
tail of the sum y',-^^ , ,^ O1O2 . . . Oj_i log — < -777- We will show how to choose niQ > nii 



to satisfy the conclusion of the lemma. 



Oi 



40 
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We bound the influence of tlie cliange from (3^ to (3^ using Lemmas 13.41 and K15I The 
influence on each of the "head elements" {i < n + rrii) is bounded by 



ai(/5i)...a,_i(/5i)log 



log 



a^m 



ai(/?^)...ai_i(/?^)log 



ai{(3 



N\ 



i-1 



^ \ ^ r)j — (n+m)+l_|_r)i— (n+m)+2 ^ r)i— (n+m)+3 ^ omi+S- 



By making m sufficiently large (i.e. by choosing a sufficiently large mo) we can ensure that 

ai(/5^)...a,:„i(/5^)log-^ 



1- 



40$(/3i) 



< 



ai{l3 



N\ 



ai(/3i)...«,_i(/3i)lo^ 



1 



<1 + 



a,iP^] 



A0^{(3^)' 



hence 



ai(/?^^)...«,_i(/?^^)lo, 



ai(/?i)...a,_i(/?i)log 
1 



ai(/5^) "^^" .•••--.V- ^'""a.(/3i; 



< 



^^^^a,(/?^)...a._,(/?^)log^^^. 
Adding the inequality for i = 1,2, . . . ,n + rrii — 1 we obtain 

n+mi — 1 n+mi — 1 

^ ai(/?^)...a,_i(/5^)log-^^^- ^ «i(/3i)...a._i(/3i)log-^ 



i=l 



1=1 



40$(/3 



- '5^ ai(/5^)...a._i(/3^)log-i-< 



< 



n+-mi — 1 



i=l 



«,(/?!) 40<l>(/3 



TS^iP'] 



e 
40' 



Hence the influence on the "head" of $ is bounded by — . 

^ 40 
To bound the influence on the "tail" we consider three kinds of terms 

ai{P^) . . . ai-i{P^) log — ^ : those for which n + rrii <i<n + m — 2,i = m + n — 1 
and i > m + n + 1 (recall that i = n + m is not in $^). 
• For n + rrii <i<n + m — 2. By Lemmas 13.41 and 13.51 

1 



a^iP')...a,.,iP')lo^ 



log- 



aiilS^] 



ai(/5^)...«,_i(/3^)log- 



1 



j-i 



^ \ ^ rfj-{n+m)+l _|_ 2*-("+"i)+2 ^ r)j-(n-|-m)-|-3 ^ 9 



Hence in this case each term can increase by a factor of at most e^. 
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• For i = n + m — 1 Note that the change decreases log so that log < 

'-^n+m— 1 'J^n+m— ll,P ) 

log L_. hence we have 

n+m— 2 

i=i 
Hence this term could increase by a factor of at most e. 

• For i > n + m + 1: Note that aj for j > n + m are not affected by the change, and the 
change decreases a„+m, so that an+m{P^) < C(n+m{P^)- Hence 

ai(/?^)...«,_i(/?^)log ^ 



ai(/3i)...a,_i(/5i)log^— «i(/?i)...a„+„(/3i) " 

"UP ) 

/rjATN ('/^Af^ n+m~l 

So in this case each term could increase by a factor of at most e^. 

We see that after the change each term of the tail could increase by a factor of e at most. 



e'e 



, hence the change in the tail 



The value of the tail remains positive in the interval 0, 

^ V 40 

is bounded by — < — . 
^ 40 40 

So the total change in $~ is bounded by 

e 'de e 
change in the "head" + change in the "tail" < 1" 1~ ~ 7- 



D 



The following Lemma follows immediately from Lemma lo? 



Lemma 3.7. For any e and for the same mo{e) as in Lemma Vj. 61 for any m > tuq and 

N, 

|<l>-(/?^)-$-(/3^+i)|< J. 
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Proof. We have 

i$-(/5^) - <i'-(/3^+^)i < i-K-i/?^) - $-(/5^)i + \<^-m - $-(/5^+^)i < I + f = |- 



n 



We will now have to take a closer look at the term ai . . . an+m-i log = $(^7) — $ (u;) 

We will need the following simple statement. 

Lemma 3.8. For any k > 1, ak-iUk < -■ 

Proof. There is an integer / > 1 such that 

1 



1 1 

afc-ittfc = 1— Oik < ■ ctfc = -. 



1 



D 



Denote $"^(0;) = ai . . . an+m-i log = $(u;) — $ (u), we are now ready to prove the 

following. 

Lemma 3.9. For sufficiently large m, for any N , 

$i(/3^+i)-<l>i(/5^)<|. 



Proof. According to Lemma [3.31 we have 



log 



«i(/?^+i)...«„+„_i(/3^+i) 



«i(/3^)...a„+^_i(/?^) 



< 



n+m— 1 

E2 



i=l 



i-{n+m) /pj ^ 

' N' 



Hence ai(/3^+i) . . . a„+^_i(/3^+i) < ai(/3^) . . . a„+^_i(/3^)ei/^, and 



log 



$l(/3^+l) < $l(/?^)el/^ «n+>n(/?^+^) ^ $l(^^)el/iV 



log 



«n+m(/5^) 



l/^l 0g(iV+ 1 + 1/0) 

log(iV+l/0) • 



Hence 



V log(A^ + 1/0) 



< 



$^(/?^)((l + 



e \ log(Ar + 1 + 1/0) 



NJ log(iV + l/0) 



1 . 
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We make the following calculations. Denote x = — - — —- ttt — , then (A^ + l/dY 

log(Ar + l/0) 

N +1 + 1/0, and 

1 

^ ^^^ A^ + 1/0 N+l/(P 

3 3 3 

A^ + 1/0 > e^/^, and so x - 1 < — — < — , thus x < 1 + — . 

' iV + 1/0 iV A^ 

By Lemma [3.81 we have 

1 / . N {n+m-2)/2 

$^(/?^)=ai(/?^)...«„+^^_i(/?^)log 775]^ < U MAT + 1/0). 



Thus 



*.,.-) -*W<*w((iH-^)!t|^±iif-i)< 



- j log(Ar + 1/0) ((1 + e/Ar)(l + 3/Ar) - 1) < (^- j log(Ar + 1/0)-. 

14 e 

Since — G o(l/ log(A^+ 1/0)), this expression can be always made less than - by choosing 

m large enough. D 

Lemmas 13 . 71 and 13 . 91 yield the following 

Lemma 3.10. For sufficiently large m, for any N, 

$(^^+1) _ ^f3^) < e. 

Proof. We use Lemmas 13.71 and 13.91 For sufficiently large m, 

$(/3^+i) - $(/5^) < $-(/?^+i) - <l>-(/5^) + $^(/?^+i) - $i(/3^) <- + - = £. 

D 
To complete the proof of Lemma f3. II we will need the following statement. 

Lemma 3.11. 

lim $(;3^) = oo. 

Proof We will prove that liuiN^co^^iP^) = oo, this suffices, since $^(/5^) < $(/3^). By 
Lemma 13. 4t 

ai(/?^)...a„+^_i(/3^) 



log 



aiip^) . . . an+m-iiP^) 



^ \ ^ 2*~("+"^)+i ^ 2 



i=l 



hence 

«i(/3^) . . . «„+^_i(/3^) > i ■ aiif]') . . . an+m-i{f3') 



e2 
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and 

^^ ' e^ log(l + l/0) ^^ ' 
The latter expression obviously goes to cxd as A^ -^ cxd. D 

We are now ready to prove Lemma fH. II 

Proof, (of Lemma l3.H) . Choose m large enough for Lemma 13.101 to hold. Increase N by 
one at a time starting with A^ = L We know that ^{13^) = $(«) < $(a) + e, and by 
Lemma [3.111 there exists an M with $(/3*'^) > $(a) + e. Let A^ be the smallest such M. 
Then $(/?^-i) < <!>(«) + e, and by Lemma ICTH 

Hence 

$(«)+£ <<l>(/3^) <<l>(a)+2e. 
Choosing /3 = /3^ completes the proof. 

n 

We will now prove Lemma [3.21 

Proof, (of Lemma l3.2j) . The ^ in the expression for ^{u) converges, hence there is an 

mi > 1 such that the tail of the sum y^.^^ , ^ aia2 ■ ■ ■ cn-i log — < -. We will show how 

to choose mo > nii to satisfy the conclusion of the lemma. 

By Lemmas 13.41 and 13. 5^ for any (3'^ and any i <n + mi we have 



ai(,8^),,,ai_i(,8^)lof ^ 



log "■(''") 



ai{u) ...ai_i(u;)log 



i-l 



^ \ ^ fjj-(n+m)+l _|_ cy,-{n-\-m)+2 



ai(^) 

^ r)j— (ra4-m)+3 ^ Qn.+mi — (ra+mo)+3 omi— mo+3 

We can choose tuq sufficiently large so that e~^ ^ ° > 1 . So that 

2(s[ijj) 

ai{(3^) . . . ai_i{(3^) log > ( 1 " t^xtt ) "i*^^) • • • "i-i(^) ^^S ^n 

for i < n + mi. 

Now, for any f3'^ we have 

n+mi — 1 ^ 

$(/3^)> 5^ ai(/5>2(/5^)...«._i(/5^)log— -^> 



«^(/5^ 



n+mi — 1 
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^-2^) (*(-)-!) >^(-)-^- 

n 

We will also need the following lemma in the proof of the Main Theorem. 

Lemma 3.12. Let u = [ai, 02, 03, . . .] and let e > be given. Then there is an N = N{e) 
such that for any n > N we have $(w„) < ^{u) + e, where u^ = [ai, 02, ... , a^, 1, 1,1,...]. 

The proof is not hard and is similar to the proof of Lemma 13.61 We present the main steps 
in the proof. 

• There is an rriQ such that the sum of the tail elements of ^{u) is small: 

00 _. 

^a,M...a._,Mlog^^<^. 

i=mo ^ ' 

• Similarly to Lemma 13. 6[ we can use Lemmas 13.41 and 13.51 to show that for sufficiently 
large vnx > mo, n > rrii implies that 



} ' < 1 + —— 



ai[uj) ...a^-i[uj)\og- 



ai{uj) 

for all i < niQ. 

• Again by Lemmas 13. 4[ 13.51 and 13.81 to show that for any i (with a special consideration 

to the case i = n), 

ai{ujn) . . . ai_i{ujn) log — - — - < max { e"^ ■ ai{uj) . . . ai^i{uj) log — -— , 2^"*/^ 

Adding these up we get for n > mi: 

( e \ """""^ 1 



y^ max ( e^ ■ axUS) . . . ai_i(uj) log — ^^, 2^"*/^ | < 
,"1 V «i(^ 



00 ^ 00 



^ \ ' / i=rno ^ ' i=mQ 

$(^) + i + pJ^ + 2^—0/2 = $(^) + £ + 2^-"^«/l 
^ ^ 4 4 ■ e^ ^ ^ 2 

We complete the proof by choosing mo large enough so that 2"^"™°/^ < e/2. 
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4. Proof of Theorem 12.221 

Recall that r{6) denotes the conformal radius of the Siegel disk associated with the 
polynomial Peiz) = z^ + e'^'^^^z, or zero, if 6 is not a Brjuno number. 

We will argue by way of contradiction, and assume that there exists a countable union 
of sets 

UZiSi = M/Z 
such that the function r{6) is continuous on each Si. 

Main Lemma 4.1. There exist 

• a sequence of initial segments li = [ai, 02, . . . , a^i], and 

• a sequence of nested intervals 

[lo,ro] D [/i,ri] D [/g.rg] D ..., 
such that the following properties are maintained: 

(1) whenever i > j we have 

(2) r-i = r(7i), where 7^ = [h, 1,1,.. .]; 

(3) for each i > I and for every (3 = [Jj, t^^+i, tNi+2, ■ ■ ■] with r(/3) G [k, r,] we have 

(4) denote ii = Vi — li. Then 

ii > and ti < £j_i/2 for all i > 1; 

(5) for any (3 = [Ii,tN,+i,tN,+2, ...], i>l, we have 

$(/5)>$(7i_i)-2-(^-i). 

Proof of the Main Lemma. We prove the Main Lemma by induction on i. For the basis of 
induction, set Iq = [1], Tq = r(7o) and /q = ''^0/2, where 70 = [1, 1, 1, • • •]• Then for i = 
conditions (l)-(5) trivially hold. 

The induction step. We now have the conditions satisfied for some i > and would like 
to extend them to i + 1. 

Set 5* = 5*2+1 and let R be the set of all possible values of the conformal radius r{9) for 
9 E S. There are two possibilities: 

Case 1: There exist ^o > and rrio G N such that for every /? G S* of the form 

P = [Jj, 1, 1, . . . , 1, . . .] we have jr^ - r(/5)| > ^o- 

mo 

In this case, select < e < min(eo,^j/2). Set 

li+i = [Ii,l, 1, . . . , IJ, k+i = ri-e, and r^+i = r^ 

mo 
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7i+i = [li, 1, 1, . . .] = 7j. We have r(7i+i) = r(7j) = r^ = rj+i, hence conditions (HI),® and 
(JH) are satisfied. 

Suppose P = [Ii+i,tN,+^+i,tN,+i+2,---] with r(/5) G [/i+i,ri+i]. Then /3 ^ S" and © is 
satisfied. By Lemma 13.21 we can choose rriQ sufficiently large in Jj+i, so that for any f3 
beginning with Jj+i, we have $(/?) > $(7i) — 2~* thus satisfying (0). 

The complementary case is the main part of the argument: 

Case 2. For every e > and m G N we can find P E S starting with /j followed by m 

ones so that 

(4.1) Ti-e < r{l3) < n 

Choose an e such that ri — 3e > k > ri — Ae. Denote 

. 'log(ri -e) -log(ri - 26) log(ri - 2£:) - log(ri - 3£:)\ 
£0 = mm , > 0. 



Theorem 12. 121 of Buff and Cheritat says that the function 

v. e^^{e) + \ogr{e) 

continuously extends to M/Z. Due to compactness of R/Z, this function is uniformly 
continuous, and there exists a (5o > such that if \x — y\ < 6o then \v{x) — v{y)\ < Eq. 

We choose m large enough, so that for any ^ = [Jj, 1, 1, . . . , 1, . . .] we have |7i — CI < "^o and 

m 

SO that Lemma f3.2l holds for rriQ = m with I = li and e = 2~*. Write 

m 

By assumption, the conformal radius r(«) is continuous on S. Hence there is a 5 > such 
that 

\r{x) — r{P)\ < e whenever |a: — /3| < 5 and x E S. 
By Lemma f3.12t there is an A^ such that for any n> N, 

Pn = [h) I) I5 • • • 5 I) '^Ni+m+li ■ ■ ■ 1 tNi+m+nj I5 I5 • • -J 
m 

satisfies 

^Pn)<HP)+eo. 
We can choose n > N large enough so that for any x whose contnued fraction expansion 
has the initial segment 

-'« ~ i-'t) 1; 1; • • • ; 1 ; iNi+m+1, iNi+m+2, • • • 5 ^TVi+m+raJ; 
m 

we have 

\x- P\ < min((5o/2,(5). 
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Start with uq = (3n = [li, 1, 1, • • •]• We have |a;o — /3| < So, and hence |i^(co'o) —v{P) \ < Eq. 
So 

logr(cjo) =v{uJo) -$(cJo) >v{f3) - Eq - ^(3) - Eq = \ogr{(3) - 2eo. 

By Lemma 13.11 we can extend /° to a longer initial segment // so that setting tui = 
[//, 1,1,.. .] we have 

$(cJo) + 2^0 < $(^i) < $(c^o) + 4eo- 

We have Ic^o — P\ < ^o/2 and \uji — P\ < 5o/2, so \uJo — uji\< 6o, and |t;(co'o) —v{uji)\ < Eq. 
Hence 

log(r(u;i)) = v{uji) - ^(cui) > v{uJo) - Eq - ^{ujq) - Aeq = log(r(u;o)) - S^o, 
and 

log(r(u;i)) = v{uji) - ^{uji) < v{uJo) + Eq - ^{ujo) - 2eo = log(r(u;o)) - ^o- 
Hence 

log(r(u;o)) - Seo < log(r(a;i)) < log(r(cjo)) - £:o- 

In the same fashion, we can extend 1} to If, and obtain U2 = [If,l,l,...]so that 

log(r(u;i)) - 5eo < log(r(a;2)) < log(r(u;i)) - £:o- 

Recall that 

log(r(^o)) > log(r(/5)) - 2eo > log(ri ~ e) - 2eo > log(ri - 2e) + G^o- 

Hence, after finitely many steps, we will obtain /f and Uk = [If , 1, 1, . . .] such that 

log(ri - 3e) + eo< \og{r{uJk)) < log(ri - 3e) + Qeq < \og{ri - 2e). 

Choose li+i = If, 7i+i = iUk, k+i = h, and r^+i = r{uJk). We have 

k+i <ri-3E < n+i <ri- 2e. 

Condition ([H) and © are satisfied by definition. Condition (JH) is satisfied because 

ii = ri - k < Ae and ii+i < ii - 2e < ii/2. 

Condition Q is satisfied, because x = [If ,...]& S only when r(x) G [r{P) — e, r(/3) + e] C 
[vi — 2e, ri + e], and by our construction [/j+i, rj+i] is disjoint from [r^ — 2e, Tj + e]. Finally, 
condition (jS)) is satisfied because Lemma (3.21 holds with e = 2~* with the m chosen in the 
initial segment /?. We have thus completed the proof of the Main Lemma f4. 11 

D 
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Passing to the Limit. The completion of the proof of Theorem l2.22l rehes on the foUowing 
lemma: 

Lemma 4.2. Denote 7 = limj^oo7i- Then the following equalities hold: 

$(7) = lim $(7j) and r{'~f) = lim r(7j). 

i— >oo i— >oo 

Proof. By the construction, the limit 7 = lim7j exists. We also know by condition (JU of 
the Main Lemma f4. II that the sequence r(7j) = Tj converges uniformly to some number r. 
By condition ^ of the Main Lemma l4. II the sequence $(7i) — 2"*-*"^-' is non-decreasing, 
and hence converges to a value ip (a priori we could have ip = 00). The sequence $(7,) 
must converge to ip as well. 

By Proposition 12.171 we have r{'j) > r > 0. On the other hand, by condition (0) of 
the Main Lemma f4.H we know that $(7) > $(7j) — 2~'^*~^^ for all i. Hence $(7) > ^. In 
particular ip < 00. 

From |B("2j we know that 

(4.2) ip + logr = lim(<l>(7i) + logr(7i)) = $(7) + logr(7). 

Along with r{'~f) > r and $(7) > ip this yields $(7) = ip, and r(7) = r, which completes 
the proof. D 

Finalizing the argument. Let 7 be the limit from the previous lemma. We claim that 
7 ^ USj. Indeed, for every i, the continued fraction Jj is an initial segment of the continued 
fraction expansion of 7 by condition (0) of the Main Lemma 14.11 By Lemma 14.21 

r(7) = limr(7i) = limrj G [li,ri]. 

Thus by condition Q of the Main Lemma [4. II we have 'j ^ Si. We have in this way arrived 
at a contradiction with USi = M/Z, which completes the proof of Theorem 12.221 

5. Concluding remarks 

Connection with the work of Buff and Cheritat. Let us outline here how the methods 
of jBClj can be applied to prove Theorem 12.221 instead of the estimates of ^ (we note 
that a newer version of the same result exists [ABC , where the arguments we quote are 
simplified). The main technical result of that paper is the following. Let a = [ao,ai, . . .] 
be a Brjuno number, and as before denote pk/qk the sequence of its continued fraction 
approximants. Let A > 1 and for each integer ra > set 

a[n] = [ao, ai, . . . , a^, A^", 1, 1,1,...]. 

Then for this particular sequence of Brjuno approximants of a, 

$(a[n]) — > $(a)+logA, 



and moreover 



limr(Q;[n]) = r{a)/A. 
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The last equality can be used to construct the "drops" in the value of the confornial radius 
of the Siegel disk needed to inductively avoid the sets Si. In this way, one obtains a sequence 
of Brjuno numbers 6i -^ 9 with confornial radii Tj = r(^j) > rj+i such that limrj = r > 0, 
and 9i is not in any of the Sj up to z-th. 

It remains to show that r{9) = r, as a priori only the inequality "<" is known. Buff 
and Cheritat demonstrate it in their context. The idea is, roughly speaking, in showing 
that the boundary of A(6'j) is well approximated by a periodic cycle of a high period. The 
perturbation 9i h-» 9i^i is then chosen sufficiently small so that the cycle does not move 
much. 

As a final remark, let us point out: 

Remark 5.1. Combining the methods of |BClj with our argument as outlined here one may 
strengthen the Main Theorem by showing that there exists a non-computable Siegel Julia 
set for which the boundary of the Siegel disk is smooth. 

Further progress. It is a natural question to ask whether the construction of non- 
computable Julia sets carried out in this paper can be replaced with a different, perhaps, 
simpler approach. Jointly with I. Binder, we have demonstrated the following in |^BB Ylj : 

Theorem 5.1 f |BBYlj ). Let R he a rational mapping of the Riemann sphere with no 
rotation domains (either Siegel disks or Herman rings). Then its Julia set is computable 
by a TM M'^ with an oracle for the coefficients of R. 

Moreover, it is shown in the same paper that the Julia set of a quadratic polynomial Jc 
with a periodic Siegel disk with confornial radius r is computable by a TM with an oracle 
for c if and only if r itself is computable by some such machine. In retrospect, therefore, 
our approach finds the only available class of examples. 

The size of the set of parameter values 9 G M/Z for which J{Pe) is uncomputable is 
rather meagre. One can show combining the results of BBY l with, for example, those 
of Petersen and Zakeri |PZj that this set has Lebesgue measure zero; and Theorem 15.11 
implies that its complement contains a dense Gs subset of M/Z. It is natural to ask if, for 
example, its Hausdorff dimension is positive, and the answer to this question is not known 
to us. It is also interesting to ask if any values of 9 in this set are computable reals (as 
there are only countably many computable reals, and our procedure clearly produces an 
uncountable set of 6''s, most of them cannot be computable). We again do not know the 
answer to this. 

On the practical side of things, to our knowledge, one has not been able to produce infor- 
mative pictures of quadratic Julia sets with Cremer orbits, although by Theorem 15.11 this 
is theoretically possible. One potential explanation is that the computational complexity 
of these sets (the amount of time it takes to decide whether to color a pixel of size 2~"' as 
a function of n) is very high. This is indeed so for the naive algorithms. In |BBY2j jointly 
with I. Binder we have constructed quadratic Julia sets whose computational complexity 
is arbitrarily high, but again all with Siegel disks. 
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A natural first step towards studying the complexity of Cremer Julia sets is to look at 
parabolics, but the first author has recently demonstrated in |Brv2j that having a parabolic 
orbit does not qualitatively change the complexity of computing a Julia set. This opens 
an entertaining possibility that some Cremer Julia sets have attainable computational 
complexity, and could be practically drawn by a clever algorithm. 
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